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Granular fronts are a common yet unexplained phenomenon emerging during the gravity driven 
free-surface flow of concentrated suspensions. They are usually believed to be the result of fluid con¬ 
vection in combination with particle size segregation. However, suspensions composed of uniformly 
sized particles also develop a granular front. Within a large rotating drum, a stationary recirculating 
avalanche is generated. The flowing material is a mixture of a visco-plastic fluid obtained from a 
kaolin-water dispersion, with spherical ceramic particles denser than the fluid. The goal is to mimic 
the composition of many common granular-fluid materials, like fresh concrete or debris flow. In 
these materials, granular and fluid phases have the natural tendency to segregate due to particle 
settling. However, through the shearing caused by the rotation of the drum, a reorganization of the 
phases is induced, leading to the formation of a granular front. By tuning the material properties 
and the drum velocity, it is possible to control this phenomenon. The setting is reproduced in a 
numerical environment, where the fluid is solved by a Lattice-Boltzmann Method, and the particles 
are explicitly represented using the Discrete Element Method. The simulations confirm the findings 
of the experiments, and provide insight into the internal mechanisms. Comparing the time-scale of 
particle settling with the one of particle recirculation, a non-dimensional number is defined, and is 
found to be effective in predicting the formation of a granular front. 

PACS numbers: 45.70.Ht, 45.70.Mg, 47.57.Gc, 47.57.E- 


I. INTRODUCTION 

Concentrated suspensions are a prominent multiscale 
material, since they exhibit a wide range of physical phe¬ 
nomena, whose relative importance varies widely at dif¬ 
ferent length scales. The study of concentrated particle- 
fluid suspensions m, is an essential research field due 
to its many applications, among which the processing 
of food and pharmaceutical materials, the construction 
industry, and the study of natural hazards are notable 
examples. Among those, phase separation is typical of 
suspensions where the length scale is sufficiently large so 
that particles are able to develop inertia independently 
of the fluid. This has dramatic and potentially undesir¬ 
able effects, since the material loses its uniformity and 
becomes therefore problematic to process. The modeling 
of phase separation is still a challenge, both numerically 
and analytically, since constitutive relationships able to 
describe the resulting non-homogeneous material have 
difficulties achieving generality. 
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Particle suspensions with grains heavier than the fluid 
vertically segregate due to sedimentation. The same sus¬ 
pension, if flowing down an incline, can also experience 
other segregation patterns. Examples are inverse and 
lateral grading as a function of the confinement condi¬ 
tions (2D or 3D flow). In many cases, a concentration 
of particles at the front of the flow is observed. This 
behavior is well documented PHS] and can be observed 
in nature [6;. A particularly important example are de¬ 
bris flows, an avalanche that originates in mountainous 
terrain composed of a mixture of water with granular 
material [7]. For still debated reasons, the front of the 
debris flows is often richer in large boulders than the 
rest of the flow. This granular front brings a signifi¬ 
cant contribution to the destructive potential a debris 
flow represents. Similar mechanisms are also observed 
in snow avalanches 0, pyroclastic flows, and rock slides. 
At smaller scales, grain-rich fronts have been observed 
on thin film particle-laden flow mm- The formation 
of a granular front is usually explained as an effect of 
inverse size segregation 0HI]: When a mixture of parti¬ 
cles flows down an incline, the shear deformations causes 
an increase in the porosity of the medium, and a number 
of void spaces are continuously created. If the particles 
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that populate the system are not of the same size, it is 
more likely that smaller particles will fall into the newly 
created voids. As a consequence, larger particles tend to 
migrate towards the surface where they are advected by 
the flow and brought towards the front. However, the 
presence of particles of different sizes is not necessary to 
generate a particle-rich front. Kaitna et al. |12| showed 
that the same pattern can be observed with particles of 
uniform size immersed in a fluid. They studied how the 
flow characteristics influence the formation of a granular 
front, showing that the phenomenon is dependent on the 
particle content, fluid viscosity, and on the shear condi¬ 
tions. This scenario however is much less studied due 
to the difficulties in performing experiments and simula¬ 
tions. Thus the mechanism that induces the formation 
of granular fronts remains unclear. 

With the use of a rotating drum IT3ll23l (see 
Figs. 0 (a), IT stationary, free-surface recirculating flows 
are obtained, whose physical properties can be continu¬ 
ously measured. The segregation pattern can be induced 
and controlled by regulating the rotational velocity of 
the drum as well as by modifying the particle concen¬ 
tration of the mixture (see Fig. [l] (d)). In order to gain 
insight into the mechanisms that govern the transition 
to a granular front, the experimental settings are repro¬ 
duced in a numerical environment. The used method is 
based on the well-established coupling between the Dis¬ 
crete Element Method (DEM) for the representation of 
the particle phase, and the Lattice-Boltzmann Method 
(LBM) for the solution of the fluid phase. 

This paper offers a desc ripti on of the experimental 
(Sec. 0 and numerical (Sec. |III| ) setting, followed by the 
presentation of the results in Sec. |IV| and their interpre¬ 
tation in Sec. El Sec. ED summarizes and comments on 
the conclusions. 



FIG. 1: Schematics of the particle fluid mechanisms 
occurring in a rotating drum test. The color scheme is 
blue for the fluid and red for the particles, (a) Rotating 
drum reference system and common particle-fluid 
distribution, (b) Particle and fluid motion patterns, (c) 
LBM and DEM discretization mesh employed in the 
numerical simulations (a two-dimensional view is 
presented for simplicity), (d) Granular front formation 
process by increasing the rotational velocity of the 
drum uji. (e) Fluid-particle and particle-particle 
interactions considered as a result of a granular flow 
through a viscous flowing fluid. 


II. EXPERIMENTAL METHODS 

The rotating drum employed throughout this study is 
located at the University of Natural Resources and Life 
Sciences (BOKU), Vienna, Austria, and described in de¬ 
tail in Refs. 

niiiiani]. Therefore, only a brief descrip¬ 
tion of the instrumentation and the configuration of the 
drum for the experiments in this study is presented. The 
rotating drum has a diameter of R = 2.46 m and a rect¬ 
angular inner section of width W = 0.45 m (see Fig. [2|. 
The cylindrical reference system used in this paper is 
presented in Fig. [2j with the angular coordinate 0 being 
0° at the 6 o’clock position and 90° and -90° at the 3 
and 9 o’clock position, respectively. The side-walls are 
composed at one side of stainless steel and at the op¬ 
posite side of acrylic glass, while the bed of the channel 
is roughened using a rubber surface with protrusions in 
a zigzag pattern of approximately 3 mm in height and 
5 mm in separation. Material losses are observed due 
to the clogging of material in the roughened bed, which 
then drops and accumulates at the inner roofs (Fig. [2^). 


The loss of material is recognized as one of the major 
limitations of the tests conducted, and is noticeable for 
high viscous flows. The drum rotates around its axis at 
a constant rotational speed uj of approximately 0.3, 0.5, 
0.7, or 0.9 rad/s. 

The material is chosen to represent a simplified nat¬ 
ural suspension, resembling to some extent the one in 
a debris flow. Particles are ceramic balls of density 
p s = 2420 kg/m 3 and diameter d = 32.6 =b 0.03 mm, 
whose high strength and stiffness enables the occurrence 
of strong collisional forces without large deformations 
or material failure. The friction angle between the ce¬ 
ramic balls and the drum roughened bed material is 
^bed = 42.5°, measured in a similar way as presented 
in Ref. [28], while the inter-particle friction angle is 
2 p s = 27.7° [29]. Particles collide with a constant co¬ 
efficient of restitution of c = 0.7 when dry, and c ^ 0 
when covered with fluid. 

The fluid phase is obtained by mixing kaolin powder 
(particle size of 2.0 pm (D$ o), and particle density of 
2600 kg/m 3 ) with water. Changing the relative propor- 
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FIG. 2 : Rotating drum at BOKU. a) Height point-laser, 
b) Normal load cell, c) Fluid pressure transducer, 
d) High-speed camera, and e) Inner-roofs (drop 
protection). Adapted from Refs. [24l [25] . 


tions of kaolin and water, mixtures of different rheolog¬ 
ical behaviors are obtained m- The parameters of the 
kaolin-water dispersion are measured using a simple co¬ 
axial cylinder rheometer (Bohlin Visco 88 ) with a gap of 
1.5 mm. The mixture used in the experiments is com¬ 
posed of kaolin and water in equal parts by mass, result¬ 
ing in a fluid with density pf = 1420 kg/m 3 and a rheo¬ 
logical flow curve as presented in Fig. [3] For simplicity in 


the numerical analysis (Sec. Ill, IV]), the resulting dis¬ 


persion is assumed to behave like a plastic fluid with a 
yield stress To, whose constitutive relation between shear 
stress r and shear rate 7 is approximated with a Bingham 
law 


T = To + MoT (!) 

where po is the plastic viscosity (see Fig. [3|. 

During the experiments, a volume (net of losses) of 
40 kg of fluid is added to the drum, together with a vari¬ 
able amount of particles (10, 20, 30, 40, 50 kg). This 
results in mixtures with a global particle concentration 
in the range 13, 23, 31, 38 and 43% respectively. A cou¬ 
ple of initial rotations are completed to uniformly cover 
the channel bed and reach steady state. 

The drum is instrumented as shown in Fig. [ 2 ] At every 
successive rotation, basal total load and basal fluid pres¬ 
sure are measured with a set of two load cells (HBM 
PW2GC3) and one piezoresistive pressure transmitter 
(Keller PR25Y), fixed normal to the bed and displaced 
180° along the circumference. The load cells are con¬ 
nected to a plate, 60 mm in diameter, covered by the 
same roughened layer of the drum’s bed. At the same 
time, the flow height is recorded by a point-laser sensor 
(Baumer OADM 20), rotating with the drum yielding a 
complete height profile every turn of the drum. All in¬ 
struments record at a sampling frequency of 1200 Hz. In 
addition, a high-speed camera (Optronis CR3000x2) with 



FIG. 3: Rheological behavior of the kaolin-water 
dispersion. The flow curve is approximated as a 
Bingham fluid (solid line). 


a 28-200 mm lens records the front of the flowing mate¬ 
rial. The camera is mounted near the center of rotation 
of the drum, with a focal distance of 1.2 m to the channel 
bed. The videos are recorded at a frame rate of 500 fps, 
with a frame size of 1696 x 1710 pixels. The images al¬ 
low to reconstruct the front location and to identify the 
position of the particles (see Sec. IV). 

The material inside the drum forms a stationary free- 
surface flow, whose front, body, and tail present distinct 
features. Along the channel axis, the front presents ve¬ 
locities both in the 0 and r directions, while in the body, 
velocities in r are negligible. At the tail the flow gives 
rise to unsteady avalanche releases ED. 


III. OUTLINE OF THE LBM-DEM METHOD 

The experimental setting described in the precedent 
chapter is recreated in a numerical environment. The 
natural choice for the simulation of the particle-fluid mix¬ 
ture is a hybrid DEM-LBM method, due to its flexibil¬ 
ity to handle the many challenges posed by the non¬ 
linear boundary conditions. A detailed presentation of 
the DEM-LBM scheme is beyond the scope of the present 
work. The readers are referred to the literature on the 
LBM in Refs. [3^, 33] and on the DEM in Ref. [34]. The 
coupling between these two solvers has been the object 
of much recent research [35] and has nowadays achieved 
a high level of maturity, with the possibility to simulate 
real systems also with non-Newtonian fluids [36;. 

The domain is discretized in a fixed, regular, cubic 
lattice, whose nodes belong to these five categories: fluid, 
gas, interface, solid walls, or solid particle. The role and 
governing equations for every node type are described 
below. 

Fluid nodes are the locations where the fluid dynam¬ 
ics is solved. This is realized through a discretization 
of the Boltzmann equation. The fluid is not represented 
by continuous fields of macroscopic variables, like in most 
classic solvers, but rather through streams of microscopic 
particles, or populations, whose collisions and subsequent 
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FIG. 4: Snapshots of the DEM-LBM simulations 
(particle content of 31%). (a) The formation of a 
granular front by an increase in the angular velocity of 
the drum. The tangential velocity vq along a 
longitudinal profile for fluid (b) and particles (c) . 


redistribution governs the dynamics. The population 
density / at every fluid node is used to reconstruct den¬ 
sity pf, pressure pf, and velocity Uf through simple sum¬ 
mations. If fi describes a population moving at speed c*, 
this translates into 

Pi = E Pi = C l' Pi > U i = E h C i/Pi > ( 2 ) 

i i 

which implies that the fluid is actually treated like a 
slightly compressible medium, and the pressure is a sec¬ 
ondary variable obtained by multiplication of pf with the 
square of the system speed of sound c 2 s . The evolution of 
fi during a unitary time step is governed by the Lattice- 
Boltzmann equation 


fi(x + Ci,t+ 1 ) = 

fi(x,t) + Qi(x,t) + Fi(x,t,g + p), (3) 


where Fi implements the effect of gravity g and of the 
fluid-particle coupling term p, while is an operator 
describing the effect of population collisions. Commonly 
expressed by the Bhatnagar-Gross-Krook linear approxi¬ 
mation E2. it relaxes the system to an equilibrium state 

neq 

Ji ? 


a = 


fp - fi 


( 4 ) 


The transition to equilibrium is described by the relax¬ 
ation time r, which is proportional to the plastic viscosity 
of the fluid po and its yield stress To [38] as 

H+yr 2 - ® 

z c s 

Gas nodes represent the space not occupied by the 
fluid, and therefore neither contain nor transfer popula¬ 
tions. 

Interface nodes represent the interface between fluid 
and gas, and are similar to fluid nodes, in the sense that 
the streaming of population happens in an identical fash¬ 
ion. However, they are granted an additional degree of 
freedom, a variable called mass m G (0,1), used to track 
the evolution of the surface. Interface nodes mutate into 
fluid nodes if m > 1 or into gas nodes if m < 0. The 
evolution of mass depends on the difference between the 
populations streaming in and out of the node, 

m(t + 1) = m(t) + y] a (/ in - / out ) (6) 

where a depends on whether the population exchange 
happened between two interface nodes or between a fluid 
and an interface node |39]. The modeling of surface ten¬ 
sion can be included in this formulation. However for 
the studied cases a set of test simulations revealed sur¬ 
face tension to have no sensible effect on the results, and 
its modeling has therefore been neglected. 

Solid wall nodes are all nodes that he inside solid 
walls, i.e. the drum cylinder and lateral walls. They 
do not contain nor transfer populations, but neverthe¬ 
less affect the fluid since they both enforce no-slip at the 
boundary and transfer momentum to the fluid when the 
drum is in motion. No-slip is enforced by requiring all 
population streams fi pointing towards a solid node to be 
reflected back. The reflected populations fy are modified 
taking into account the momentum transfer, 

fi' = fi- 6wiPfU w ■ a, (7) 

where u w is the velocity of the wall at the reflection lo¬ 
cation and Wi is a weight depending on the lattice config¬ 
uration. Nodes belonging to the drum surface employ a 
modified version of this rule, to take into account the cur¬ 
vature of the boundary. Details about the curved-surface 
treatment can be found in Refs. mm- 


TABLE I: Material parameters used in the DEM. 


Parameter 

Value 

Mass density 

2420 kg/m 3 

Young’s modulus 

1.6 • 10 7 Pa 

Poisson’s ratio 

0.3 

Restitution coefficient 

~ 0 

Shear damping 

10 kg/s 

Particle-particle friction angle 

27.7° 

Particle-drum friction angle 

42.5° 
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Experiment 




FIG. 5: Experimental and numerical flow height (5a) 
and bulk basal pressure profiles (5b) at 0.3 rad/s and 
particle concentration of 31%. The dashed line has the 
height of one particle diameter. 


Solid particle nodes are all nodes contained inside 
solid particles. For the particle-fluid coupling algorithm 
the direct forcing approach m is employed in a sim¬ 
plified form similar to the one described in Ref. E5) . A 
sketch of the interaction mechanism is given in Fig. [11(c). 
Particles are immersed in the fluid and are advected 
through the LBM regular grid. The difference between 
particle and fluid velocity is used to compute a hydro- 
dynamic interaction force that is transmitted both to 
the fluid and to the particle equations of motion. The 
Lattice-Boltzmann equation is solved in this case with 
the additional forcing term p appearing in Eq. [3J calcu¬ 
lated as 


p(x, t) = Pf(x, t) [u{(x, t) - u p (x, t )], ( 8 ) 

where u v is the velocity of the particle at the node lo¬ 
cation [44]. A force of opposite sign is applied to the 
particle. 

Particle dynamics is solved using a DEM model (see 
Fig. [T] (e)). The collisions between particle, and the con¬ 
sequent momentum exchange is modeled using a Hertzian 
law [43j 05] with constant coefficient of restitution. The 
material parameters used in the simulations are recapped 
in Table IT] 


IV. COMPARISON BETWEEN NUMERICAL 
AND EXPERIMENTAL RESULTS 

The height and the basal pressure longitudinal profiles 
obtained in experiments described in Se c. [TT| are used to 
validate the DEM-LBM scheme. Fig. |5a| presents the 
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FIG. 6: Front as a function of the angular velocity of 
the drum and particle concentration. 


height profile as detected by the point-laser sensor and its 
numerical counterpart. Both profiles are averaged over 
all available data in steady state. The experimental pro¬ 
file is based on 6 recorded data sets, one for every full 
rotation of the drum, while the numerical profiles are ob¬ 
tained from a much larger set of data. Good agreement 
is found between the numerical simulation and the ex¬ 
perimental measurements, matching reference points as 
the front location and flow height. In a similar way, the 
experimental basal pressure profile (Fig. |5b| ) coincides 
with the numerical pressure, accounting simultaneously 
for the particle and fluid contributions. 

The position of the particles inside the fluid cannot 
be clearly obtained with the sensors. Therefore the se¬ 
quence of frames obtained with the high-speed camera 
(see Fig. [6]) is analyzed with an edge-detection tech¬ 
nique m for the flow front recognition. The most strik¬ 
ing characteristic of the flow is the possibility to obtain a 
granular front, even with a relatively low initial concen¬ 
tration of the particles. The higher the rotational speed 
of the drum, the more particles will get close to the fluid 
front. The formation of a granular front strikingly re¬ 
semble the one observed in a debris flow. The process is 
shown in Fig. [7] where the fluid front position 0 f, as iden¬ 
tified by the sensor is shown with blue markers, while the 
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FIG. 7: Phase front transition. Experimental and numerical simulation results are shown. The experimental fluid 
front is extracted from the recorded height profiles, while the experimental particle front comes from the analysis of 

the the images obtained with the high-speed camera. 


particle front position 0 P as obtained by the high-speed 
camera, is shown with red markers. Fig. [7] also shows 
the position of the center of volume (CoV), obtained by 
integration of the height profiles. The values obtained 
from the simulations are shown with continuous lines. 

In agreement with Ref. [ 12 ], the transition from a fluid 
front into a particle front can be induced by an increase in 
the angular velocity of the drum (see Fig. [ 6 |. At low an¬ 
gular velocities, the particles are concentrated at higher 
inclinations. Then, with increasing velocities at the base, 
the particles diffuse, and start to move towards lower in¬ 
clinations. Along with the motion of the particles, the 
interstitial fluid is strongly sheared by the particle mo¬ 
tion and drum bed. As a result of this shearing, the fluid 
front moves backwards while the particles move forward 
until the particle front is formed. 


V. INTERNAL MECHANISMS GOVERNING 
PHASE SEPARATION 

The emergence of granular fronts, described in the pre¬ 
vious section, is the result of two mechanisms. The first is 
the migration of the fluid front towards higher positions, 
observed for growing angular velocity of the drum, and 
the second is the movement of the particles under high 


shear. Simple analytical models for these two phenomena 
are described in the following. 

As seen in Fig [3] the Kaolin-water mixture can be 
approximated with a Bingham fluid, characterized by a 
yield stress To and a plastic viscosity /ro- Obtaining an 
analytical solution of the Navier-Stokes equations for a 
Bingham plastic in the partially-filled rotating drum is 
not an easy task: The fluid presents a free surface whose 
shape is not given a priori. However, considering a single 
slice of the drum (6 fixed) and approximating it with a 
one dimensional flow in the 0 direction, an approximate 
solution can be obtained, which is analogous to the use 
of a long-wave approximation m ■ With this hypothe¬ 
sis the flow is assumed to be steady and self-similar over 
the width w of the drum and the edge effects due to the 
side walls of the drum are neglected. The shear rate and 
shear stress tensors each boil down to one component 
7 Qr = 7 ^? and tq v . The Bingham constitutive law can be 
written in this new system as: 

17^|=0 if (\t9v\ < t 0 ), 

T 6r = (plj + Mo^) 7^ if {\rer\ > To). ^ 

The integration of the momentum equation with limits 
r = R at the base of the drum, and r = 77 at the free sur- 
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face position, gives the solution of the tangential-velocity 
profile ve(r) 

ujR — a sin 6{r — R)(r + R — 2r 0 ) if r 0 < r < R, 
ujR — asm6(R — ro) 2 if < r < ro, 

( 10 ) 

where cj is the rotational speed of the drum, and a = 
Pf^/2/io- Between the positions ro and a plug flow 
forms, whose total width ho = ro — r^ is proportional to 
the yield stress as: 


h 0 = 


rp 

Pig sin 0 ' 


(ii) 


By imposing stationary conditions, and therefore total 
flux qo = 0 in the 0 direction, the total height of the fluid 
h = R — rh is obtained. The following expression can be 
obtained by integrating the velocity profile (Eq. [To| ) over 
the whole height h (both yielded and unyielded regions): 


1 9 

qe = uoRh - asm6(h — ho) (2h + ho) = 0. (12) 

O 


This equation can be solved numerically to obtain the 
height h of the flow. An elegant solution can be extracted 
with the approximation 2h + ho ~ 3ft, which is justified 
by ho being of the same order of h. Thus 


h — ho + 


u:R 
a sin# 


1/2 


(13) 


This solution does not hold for positions in the flow 
where a non-negligible r component exists, like the flow 
front. In addition, it is a general solution that does not 
take into account the actual flowing volume. However, 
for positions in the drum sufficiently far from the front 
(0 > 10°), Eq. [13] correctly captures the trend of variation 
of the height as linear function of the square root of the 
drum rotational speed This is so because at higher 
basal velocities, in order to keep a stationary motion, the 
fluid needs to accumulate mass by increasing its height, 
an effect typical for free-surface flows of shear-dependent 
fluids. Fig [5] shows the fluid height meas ured in e xperi- 
ments and simulations as a function of y/uj/ cj max (with 
Wmax = 1.0 rad/s). The average height h in the range of 
25° < 0 < 30° is plotted, showing a good agreement with 
the linear trend predicted by Eq. [13] 

Since the flowing volume is constant, the rise of the 
fluid height given by Eq. [13] must be accompanied by a 
transfer of mass from the front. This explains why the 
fluid front moves to higher positions when the angular ve¬ 
locity of the drum is growing (see Fig. [ 7 ]). The increase 
in the front position #f, follows the same proportionality 
to the square root of the angular velocity of the drum 
yji J, see Fig. [ 9 J as long as the influence of particles is 
minimal. In a fluid-only scenario (black solid markers), 
the relationship #f oc y/Zj holds true for the whole range 
of simulated velocities. For simulations with particles, 
there is an initial linear trend, followed by an abrupt 


0.06 

0.05 

0.04 

0.03 

0.02 

0.01 

0 

0 


FIG. 8: Analytical flow height of a free surface flow 
down an incline against the experimental and numerical 
height (average height between 25 to 30°). 
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variation. This corresponds to the point when the par¬ 
ticles get closer to the fluid front, which then becomes 
dominated by particle dynamics. For the limit of a very 
high particle content the linear trend is completely ab¬ 
sent, and the behavior of the front can only be explained 
through the behavior of particles. 

At growing angular velocities of the drum the particle 
front moves in opposite direction to the fluid fronts, to¬ 
ward lower positions. This counter-intuitive effect can be 
explained by looking at the behavior of particles located 
close to the front. The motion of these particles is gov¬ 
erned by the relative magnitude of two time-scales. The 
first is the time of settling of a particle from its suspended 
state to the bottom layer, £ S ettimg, and the second is the 
time it takes for the same particle to travel the distance 
between the grain bulk and the fluid front, travel- To 
get an estimate of Settling, the solution of the equation of 
motion for an immersed spherical particle subjected to 
gravity and buoyancy [481450] is considered: 

7T od Ujj' 7T . * 7 o 7T /-t a \ 

g Ppd = g (pp - p{)gd - —d/JoUr- (14) 

Here the flow is assumed to have a particle Reynolds 
number sufficiently low as to allow to use Stokes’ law 
for the drag force, and the fluid plastic viscosity po is 
used. The coefficient a takes into account the presence 
of surrounding particles. Its precise assessment is not 
easy in the drum condition. However, considering the 
particle packing to be close to the maximum, a can be 


estimated to be of order 10 


-2 


The solution of Eq. 14 
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FIG. 9: Numerical location of the fluid front 0f as a 
function of the angular velocity of the drum uj. Markers 
represent different particle contents. For a fluid-only 
simulation the trend is linear (solid line). 


FIG. 10: Angle between the particle bulk and the fluid 
front A 6 as a function of the time-scale ratio £, 
following a lineal trend (continuous line). Simulations 
where a granular front is observed (A 6 < 4°) are 
indicated by white markers. 


returns a settling velocity under stationary condition of 

(15) 

3 /To 

A particle has to fall through a distance equal to the 
particle layer height , which can be assumed to be pro- 
portinal to \JVpjW, with V p the total volume of parti¬ 
cles and W the drum width. The total settling time can 
be therefore estimated as 


h p 3 Mo VVp/W 

^settling — Ur — 2 a(p p — pi)gd 2 ' 


(16) 


The same particle has a velocity in the 6 direction of or¬ 
der uq = ujR (see Fig. Qb)), which is enforced by the 
boundary conditions and can therefore be considered to 
be independent from the viscosity of the fluid. This par¬ 
ticle can potentially travel the distance between the par¬ 
ticle bulk and the fluid front, Ax. This is related to the 
amount of free fluid, i.e. the fluid that is not trapped 
in the particle pores, AG = Vf — (1 — </>)Vp, where </> 
is the particle volume fraction, which can be assumed 
to be cl ose to i ts maximum value, </> max = 0.55. Since 
Ax oc yj AV/W , the travel time can be estimated as 


_ Ax _ y/AVjW 

Gravel — V, 

Uq UjR 


(17) 


The ratio between these two time-scales gives a dimen¬ 
sionless quantity able to express the tendency of a flow 


to develop a granular front. The definitions above yield 


£ = 


Gravel 

Gettling 


I AV a(p v - pi)gd 2 


V n 


fio ujR 


(18) 


When Gettling is bigger than Gravel, and therefore £ is 
small, particles will have time to reach the fluid front fol¬ 
lowing the fluid streamlines, before the settling process 
leads them to the bottom of the flow and back into the 


particle bulk. Fig. 10 shows the angle between particle 
front and fluid front A 6 = 0f — 0 P plotted as a function 
of £. On this plot, the numerical data points aggregate 
around a line, which shows that the definition of £ well 
captures the trend leading to the development of a granu¬ 
lar front. Simulations with granular fronts (A 6 < 4°) are 
indicated by white markers, and correspond to a time- 
scale ratio lower than a threshold £ t h which lies in the 
range from 0.2 to 0.4 for this specific geometry. 


VI. CONCLUSIONS 

We presented a series of experiments to study the be¬ 
havior of a mixture of fluid and particles inside a rotating 
drum. Particles are heavier than the fluid and therefore 
have the natural tendency to settle. This is contrasted 
by enforcing a recirculation of the material through the 
rotation of the drum. A sufficiently high angular veloci¬ 
ties of the drum produces a front dominated by particles. 
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This reorganization pattern of fluid and particle phases 
is typical for suspensions, both natural (e.g. debris flow) 
and industrial (e.g. fresh concrete). The same results are 
obtained in a numerical setting, using a hybrid method 
that combines LBM for the solution of the fluid phase 
with DEM for particle motion. Experiments and simula¬ 
tions show excellent agreement. 

The segregation pattern is dependent on two phenom¬ 
ena. The first is the migration of the fluid front towards 
higher positions, which can be observed for growing angu¬ 
lar velocities of the drum. This is explained by looking 
at the analytical solution for the free-surface flow of a 
Bingham fluid down an incline. At the same time, the 
particle front moves toward lower positions. The ratio 
between the time-scale of settling and the time-scale of 
particle recirculation defines a dimensionless quantity, £, 
that controls the formation of a granular front. The dis¬ 
tance between particles and the fluid front is governed by 
£, which is proportional to the square root of the global 
particle content, and inversely proportional to the an¬ 
gular velocity of the drum. This result is confirmed by 
the measurements. When plotting £ against the distance 
between particles and fluid front, numerical and experi¬ 
mental data collapse on a line. The time-scale ratio £ is 
proportional to the square of particle size and inversely 


proportional to the viscosity of the fluid, a tendency al¬ 
ready observed in Ref. m- We were not able to test 
these further dependencies, and leave the proof to a later 
work. The validity of the £ scaling can be pivotal for 
the understanding of the the free-surface flow of particle 
suspensions. In particular, it can lead to better under¬ 
standing of the condition under which a debris flow can 
develop, or a fresh concrete mixture will segregate during 
casting. 
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